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Abstract 

Considering a Vaidya exterior spacetime, we study dynamical models of prototype gravastars, 
made of an infinitely thin spherical shell of a perfect fluid with the equation of state p = a, enclosing 
an interior de Sitter spacetime. We show explicitly that the final output can be a black hole, an 
unstable gravastar, a stable gravastar or a "bounded excursion" gravastar, depending on how the 
mass of the shell evolves in time, the cosmological constant and the initial position of the dynamical 
shell. This work presents, for the first time in the literature, a gravastar that emits radiation. 
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I. INTRODUCTION 



Gravastar was proposed as an alternative to black holes. The initial model of Mazur and 
Mottola (MM) consists of five layers: an internal core < r < r 1; described by the 
de Sitter universe, an intermediate thin layer of stiff fluid r\ < r < r 2 , an external region 
r > r 2 , described by the Schwarzschild solution, and two infinitely thin shells, appearing, 
respectively, on the hypersurfaces r — r\ and r = r 2 . The intermediate layer is constructed in 
such way that T\ is inner than the de Sitter horizon, while r 2 is outer than the Schwarzschild 
horizon, eliminating the existence of any horizon. Configurations with a de Sitter interior 
lave long history which we can find, for example, in the work of Dymnikova and Galaktionov 
2|. After this work, Visser and Wiltshire [3j (VW) pointed out that there are two different 
types of stable gravastars which are stable gravastars and "bounded excursion" gravastars. 
In the spherically symmetric case, the motion of the surface of the gravastar can be written 



in the form 



3], 



he + v(r) = 0, (i) 



where R denotes the radius of the star, and R = dR/dr, with r being the proper time of 
the surface. Depending on the properties of the potential V(R), the two kinds of gravastars 
are defined as follows. 

Stable gravastars: In this case, there must exist a radius ao such that 

V(R ) = 0, V'(R ) = 0, V"(Ro)>0, (2) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such a radius Ro for which the above conditions are satisfied, the 
model is said to be stable. Among other things, VW found that there are many equations 
of state for which the gravastar configurations are stable, while others are not js]. Carter 
studied the same problem and found new equations of state for which the gravastars are 

n n n 

stable 4| , while De Benedictis et al [5| and Chirenti and Rezzolla |6| investigated the stability 
of the original model of Mazur and Mottola against axial-perturbations, and found that 
gravastars are stable to these perturbations, too. Chirenti and Rezzolla also showed that 
their quasi-normal modes differ from those of black holes with the same mass, and thus can 
be used to discern a gravastar from a black hole. 
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"Bounded excursion" gravastars: As VW noticed, there is a less stringent notion of 
stability, the so-called "bounded excursion" models, in which there exist two radii a± and a 2 
such that 

V(R!) = 0, V'ORi)<0, V(R 2 )=0, V'(R 2 )>0, (3) 
with V(R) < for R £ (R u R 2 ), where R 2 > R u 



Lately, we studied both types of gravastars |7|-[12|], and found that, such configurations 
can indeed be constructed, although the region for the formation of them in the phase space 
is very small in comparison to that of black holes. 

Based on the discussions about the gravastar picture some authors have proposed alter- 



native models 



13]- 18] . In addition, since in the study of the evolution of gravastar there is 



a possibility of black hole formation, we can find some works considering the hypothesis of 
dark energy black hole Q Q Q H H • 

In the last four years we have adopted a different approach (from VW [3]), which means 
that we started from an equation of state and found the potential of the shell jlljs]. Gen- 
eralizing the exterior spacetime in order to include a cosmological constant, we study the 
gravastar model in a de Sitter-Schwarzschild spacetime, which allowed to investigate the role 
of cosmological constant in its evolution 
a de Sitter-Reissner-Nordstrom exterior spacetime 
effects of changing the interior of the gravastar, fil 



10]. Following this direction, we also considered 



111 ] . On the other hand, we studied the 



ingit with anisotropic fluids, which can 
9fl[12|- 



be characterized by different kinds of dark energy 

Here we are interested in the study of a gravastar model whose interior consists of a de 
Sitter spacetime and an exterior radiative Vaidya's spacetime. The paper is organized as 
follows: In Sec. II we present the metrics of the interior and exterior spacetimes, and write 
down the motion of the thin shell in the form of equation ([1]). In Sec. Ill we study the model 
by using the small radiating source approximation. In Sec. IV we discuss the formation of 
black holes and gravastars, when the mass of the thin shell increases, while in Sec. V we 
study the case where the mass of the thin shell decreases. Finally, in Sec. VI we present our 
main conclusions. 
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II. DYNAMICAL THREE-LAYER PROTOTYPE GRAVASTARS 



The interior spacetime is described by the de Sitter's metric given by 



dsf = -fd? + f l dr 2 + r 2 dVt\ 



(4) 



where / = 1 - (r/L) 2 , L = ^/3/A and dtt 2 = d6 2 + sm 2 {6)d(j) 2 . 
The exterior spacetime is given by the Vaidya's metric 



ds z e = -Fdv 2 - 2drdv + r 2 dQ 2 , 



(5) 



where F — 1 — 2rn ^ , The metric of the hypersurface on the shell is given by 

dsl = -dr 2 + R 2 {r)dQ 2 , 

where r is the proper time. 

Since ds 2 = ds 2 = we find that r s = r E = R, and 



(6) 



1 d2 



ff - r'R 



F + ™ 

V 



1, 
1, 



(7) 
(8) 



where the dot denotes the ordinary differentiation with respect to the proper time. On the 
other hand, the interior and exterior normal vectors to the thin shell are given by 



n 



;-iU o,o), 



n e a = (-R,v,0,0). 



(9) 



Then, the interior and exterior extrinsic curvatures are given by 



[(3L 4 i? 2 - LH 2 + 2L 2 R 2 t 2 - RH 2 )Rt - (L + R)(L - R)(Rt - Rt)L A ] x 



(L + R)~\L - R)- 1 L 



1 r-4 



K l 0g = i(L + R)(L - R)L- 2 R 
A^ = A^sin 2 (#), 
K e TT = v 2 (2m 2 v - 3mRR - mRv + mR 2 v)R~ 3 , 
K e de = -v(2m -R) + RR, 



(10) 

(11) 
(12) 

(13) 
(14) 



Since 



23| 



[Kee] = K e de - K\ Q = -M, 



where M is the mass of the shell, we find that 



M = v(2m -R) + i(l - 2aR 2 )R. 
Then, substituting equations (0) and (JS]) into (TTT)) we get 



(15) 



(16) 



(17) 



M = -RR + {2m-R)v + R 



ir + i 



1/2 



(18) 



In order to keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a fluid with the equation of state, p = a, where a and p denote, respectively, the surface 
energy density and pressure of the shell. Then, the equation of motion of the shell is given 



by 



23] 



M + 8ttRR P = 47rR 2 [T a pu a nP} = AttR 2 fe<7ij - T ( 



a/3 U ? n i I ■ 



(19) 



where u a is the four-velocity. 

The exterior energy-momentum tensor is given by 



where 



T a(3 - 



7 £V 7 lOL f\ 

L a , L a L U, 



) 2 dm 



AnR 2 dv 



Since 



(in this paper we shall choose the plus signal), we find 
Since the interior spacetime is de Sitter, we get 



(20) 

(21) 
(22) 

(23) 
(24) 



M + 8nRRa = -mi) 6 . 



(25) 
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In order to solve this equation, let us assume that 

mi) 3 = k = kiM, 



(26) 



where k$ and k\ are constant. Note that ko < since e > 0. Thus, recalling that a 
M/(AtcR 2 ), we find that the solution of equation (1251) is given by, 



2 



M = kR~~ } (27) 

where k is a positive integration constant. Note also that k\ — — 1 is not allowed. If k\ <C 1 
we get the approximation of small emission of radiation. If k\ < then from equation (f26|) 
we can see that M > 0, meaning that the mass of the shell is increasing. In this case, 
the interior de Sitter energy transfers to the thin shell, feeding the outgoing radiation. If 
k\ > then from equation (|26|) we can see that M < 0, meaning that the mass of the shell 
is decreasing. In this case, the thin shell looses mass in order to maintain the outgoing 
radiation. Thus, for three particular values of kx, which cover all the possibilities described 
above, we have 

kR~ 2 , if taking the limit ki — > 0, 
M = I kR+ 2 , if k x = -2, (28) 
kR-\ ]lh = +1, 
which will be studied in more details in the next three sections. 

Substituting equation (1271) into equation (1T81) . we find that the potential in the form of 
equation (CQ) can be written as 

V{R, m, L, k) = — (-4:R 5 mL 2 - AL i R 2 k 2 R~^ I + 

8L*R 2 k 2 R~^ V 

2R 4 k 2 R~^L 2 + AR 2 m 2 L A + ARmL A k 2 R~^ + R 8 + k A R~^L 4 ^) . (29) 
Rescaling m, L and R as, 

fci+i 

m — > mk k i +:i J 

L ->■ Lk k i+ 3 , 

fci+i 



R ->• Rk~ A , (30) 
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we find that equation ( 129]) can be written as 

1 



V(R,m,L) 



8L 4 i? 2 



/ 9+5*:! Q 3+fcl ,3+2ki 4 

l-AtV^mL 2 - AL 4 R 2 + 2R 4 L 2 + AR 2 ~^m 2 L 4 + ARmL 4 + iT + R~~^L 4 

(31) 

Clearly, for any given constants m and L, equation ( 13TT) uniquely determines the collapse of 
the prototype gravastar. Depending on the initial value Rq, the collapse can form either a 
black hole, a gravastar or a de Sitter spacetime. In the last case, the thin shell first collapses 
to a finite non-zero minimal radius and then expands to infinity. To guarantee that initially 
the spacetime does not have any kind of horizons, cosmological or event, we must restrict 
Ro to the range, 

2m<R < L, (32) 

where R is the initial collapse radius. 

III. GRAVASTARS/BLACK HOLES WITH SMALL EMISSION OF RADIATION 

Here we will consider the limit k\ — > 0. Then, we find that 
V(R, m, L) = l — - (-AL 4 R 6 + 2R 8 L 2 + R m - AR 13 mL 2 + AR 10 m 2 L 4 + AR 5 mL 4 + L 4 ) , 



(33) 



fro which we find that lim^o V(R, m, L) = liniR_ s . 0O V(R, m, L) = —00. 
The first derivative of the potential, equation (133|) . is given by 

dV 1 



dR ~~ AL"W 



(2R 8 L 2 + 5R W - UR 13 mL 2 + 8R 10 m 2 L 4 - 2R 5 mL 4 - 3L 4 ) . (34) 



Thus the solutions for — Q are 



mi = 8R*L 



^ (7R 8 + L 2 + V9R 16 - 2R*L 2 + 25L 4 ) , (35) 



and 



m- 2 -^-^ (7R 8 + L 2 - V9R 1& - 2R 8 L 2 + 25L 4 ) . (36) 



8R 5 L 

Note that m-i and m\ are always positive, as can be seen from Fig. [TJ 
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FIG. 1: Case k\ — 0. The masses mi (left) and 1712 (right) where the first derivative of the potential 
V(R) is zero. We note that both masses are positive. 




FIG. 2: Case k\ ~ 0. The second derivative of the potential 4^-(i?,m, L) calculated at m = m\ 
(left) and m = rri2 (right). We note that ^jp(R,m = m\,L) can be positive or negative (the 
frontier between the two regions is given by equation ([39]) ) and that ^^(R,m = rri2,L) is always 
negative. 
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R 



(b) 



FIG. 3: The possible type of potentials. 

The second derivative of the potential is given by 
(hV :(R, m, L) = -t^tj (2R 8 L 2 + 45i? 16 - 8AR 13 mL 2 + 2AR w m 2 L 4 + AR 5 mL 4 + 21L 4 ) . 



dR 2 



AR 8 L 4 



(37) 



Substituting equation (|35|) into equation fl37|) we have 
;{K, m = mi, Li) — . _ . _» x 



dR 2X ' ' ' IQL A R 8 

(-2R 8 L 2 - 27R 16 - 21RW9R 16 - 2R 8 L 2 + 25L 4 + 125L 4 + 5L 2 V9R 16 - 2R 8 L 2 + 25L 4 ) . 

(38) 

Solving jj£{R, m — mi, L) = we get 

L) « 0.9731857906.R 8 . (39) 
Substituting equation ( |36|) into equation ( )37|) we have 
- — -(K, m = mo, L) = — - x 

(-2R 8 L 2 - 27R m + 21R 8 V9R 16 - 2R 8 L 2 + 25L 4 + 125L 4 - 5L 2 V9R 16 - 2R 8 L 2 + 25L 4 ) . 

(40) 

Thus, we can see from Fig. [2] that the second derivative of the potential is always positive 
at m = rrii and negative at m — m 2 . This means that the form of the potential is given by 
Figs. |3k and [3b. 
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FIG. 4: Case k\ — 0. The potential V(R,m, L) calculated at m = m\ (left) and at m = m,2 (right). 
We note that both potentials are negative. 

Substituting equation ( 1351) into equation ( 133]) we have 
V(R, m = mi, L) = -; — - x 

(AQR S L 2 + 25L 4 + 5LV9i? 16 - 2i? 8 L 2 + 25L 4 - 32 J R 6 L 4 + 9i? 16 + 3R 8 V9R m - 2R 8 L 2 + 25L 4 

(41) 

and substituting equation (1361) into equation (1331) we get 

^( J R,m = m 2 ,L) = -^^x 

^46i? 8 L 2 + 25L 4 - 5L 2 V9R W - 2R 8 L 2 + 25L 4 - 32 J R 6 L 4 + 9.R 16 - 3R S \ / 9R 16 - 2R 8 L 2 + 25L 4 ) . 

(42) 

Thus, from Fig. H]we can see that the potential is always negative at both m = mi and 
m = rn,2. This means that the form of the potential is given by the lowest curves of Figs, 
and Eb, respectively. Hence, the collapse can either form gravastars or black holes. 



IV. GRAVASTARS/BLACK HOLES WHEN THE THIN SHELL MASS IN- 
CREASES 

Now, let us assume that k\ = —2. Then, we find that 



10 



FIG. 5: Case k\ ~ 0. The potential V(R,m, L) calculated at m = mi, R c = 2 and L c = 5. This 
represents the formation of a "bounded excursion" gravastar. 



V(R, m, L) = -J— (-4i? 4 L 4 + 2L 2 R e + 4m 2 L 4 - AmL 2 R 3 + 4R 3 L*m + R 6 + R 6 L*) . 

SL R 

(43) 

Note again that \im R ^ V(R, m , L) = lim^oo V(R, m, L) = — oo. 
The first derivative of the potential, equation ( 1431) . is given by 
dV 1 



dR AL 4 R 5 



(2L 2 R 6 + 2mL 2 R 3 - 2R 3 L 4 m + R 6 + R 6 L 4 - 8m 2 L 4 ) . (44) 



Thus, the solutions for ^ = are 



R 3 

and 



R 3 



m 2 



(l-L 2 + V9 + 14L 2 + 9L 4 ) (45) 

2 (l - L 2 - V9 + 14L 2 + 9L 4 ) . (46) 

Note from Fig. [8] that m2 is always negative, while mi is always positive. Since the mass 
must be always positive, thus the unique reasonable solution for ^ = is given by mi. 
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FIG. 6: Case k\ — 0. The potential V(R,m, L) calculated at m = mi, R c = 2 and L c = 
5.324211511. This represents the formation of a stable gravastar. 

The second derivative of the potential is given by 

^(R, m, L) = -j^ji (2L 2 R 6 - AmL 2 R 3 + AR 3 L A m + R & + i? 6 L 4 + 40m 2 L 4 ) . (47) 

Thus, from Fig. |9]we can see that the second derivative of the potential is always negative 
at m = mi. This means that the form of the potential is given by Fig. [3b. 
Substituting equation (|4"5|) into equation (j47|) we have 

dR^ m = m ^ L) = -mA x 

(l4L 2 + 9 + V9 + UL 2 + 9L 4 + 9L 4 - L 2 V9 + 14L 2 + 9L 4 ) . 

(48) 

Substituting equation ( l4"5j) into equation ( 14U]) we have 
V{R,m = m u L) = x 

(-32L 4 + 30R 2 L 2 + 9R 2 - 3R 2 V9 + UL 2 + 9L 4 + 9L 4 R 2 + 3R 2 L 2 V9 + UL 2 + 9L 4 ) . 

(49) 
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FIG. 7: Case k\ ~ 0. The potential V(R,m, L) calculated at m = mi, R c = 2 and L c = 10. This 
represents the formation of a black hole. 

We notice that V(R, m = mi, L) can be positive or negative, depending on the radius R 
and the cosmological constant L (see figure ITU]) . This means that there exist possibilities of 
formation of both gravastars and black holes. 



A. Total Gravitational Mass 



In order to study the gravitational effect generated by the two components of the gravas- 
tar, i.e., the interior de Sitter and the thin shell in the exterior region, we need to calculate 
the total gravitational mass of a spherical symmetric system. Some alternative definitions 



are given by 
is given by 



24| , 25j and 26| . Here we consider the Tolman's formula for the mass, which 



Mi 



a 



Ro r n 

7T ^0 



2t 



T%drd6d<j), 



(50) 



13 



FIG. 8: Case k\ = —2. The masses mi (left) and 7712 (right) where the first derivative of the 
potential V(R) is zero. We can note that m\ is always positive and 777,2 is always negative. Thus, 
we have only m\ as solution of 4j£ = 




FIG. 9: Case k\ = —2. The second derivative of the potential ^p(R,m, L) calculated at m = m\. 
We note that ^^(R,m = mi,L) is always negative. 
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FIG. 10: Case k± = —2. The potential V(R,m, L) calculated at m = mi. We note that V(R, 
mi,L) can be positive or negative. 



1 1.5 2 2.5 3 3.5 4 
i_i i i_i i i_i i i_i i i_i i i_i i i_i i i_i i i_i i i_i i 




FIG. 11: Case k\ = —2. The potential V(R,m, L) calculated at m = mi, R c = 2 and L c = 
This represents the formation of a black hole. 
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where \J —g is the determinant of the metric. For the special case of a thin shell we have 

M G = / / ^T^5{v-R Q )dvd6d(j>. (51) 

Thus, the Tolman's gravitational mass of the thin shell is given by 

M^ heU = 3M, (52) 

and for the interior de Sitter (dS) spacetime we have 

M d G s = ~A<Rl (53) 

Thus, the de Sitter interior presents a negative gravitational mass, since Aj > 0, in agreement 
with its repulsive effect. 

Now we can write the total Tolman's gravitational mass of the gravastar as 

M total = M shell + M dS =3M _ Ia^. (54) 

This mass should also represent the Vaidya exterior mass (m = Mq 10,1 ) of the gravastar. 
This last equation can explain how the mass of the shell can increase with the time. Since 
m must decrease with the time because of the emission of radiation, the unique way that 
M may increase with the time is that the radius Rq is increasing with the time. 



V. GRAVASTARS/BLACK HOLES WHEN THE THIN SHELL MASS DE- 
CREASES 

Now, let us assume that k\ = +1. Then, we find 



V(R, m, L) = —j (-4L 4 i? 4 + 2R 6 L 2 + R 12 - 4R 9 mL 2 + 4R 6 m 2 L 4 + 4L 4 mR 3 + L 4 ) . 

oLi R 

(55) 

Note again that lim^o V(R, m, L) = lim^oo V(R, m, L) = — oo. 
The first derivative of the potential, equation (|55l) . is given by 
dV 1 



dR 2L 4 R 5 



(R 6 L 2 + 2R 12 - 5R 9 mL 2 + 2R & m 2 L A - L 4 mR 3 - L 4 ) . (56) 



Thus the solutions for ^ = are 



mi = AR*L 



l — (hR Q + L 2 + V9R™ + 2W>V + 9L 4 ) , (57) 
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FIG. 12: Case ki = +1. The mass m.2, where the first derivative of the potential V(R) is zero. We 
note from equation (|571) that mi is always positive. However, mi can be positive or negative. 



and 



1 



5R 6 + L 2 - V9R 12 + 2R e L 2 + 9L 4 



(58) 



m 2 = 



4R 3 L 2 



We note from equation (1571) that m\ is always positive and from Fig. [12] that m 2 may be 
positive or negative, depending on the radius R and the cosmological constant L. 
The second derivative of the potential is given by 



Thus, we can see from Fig. [T3] that the second derivative of the potential can be positive 
or negative at m — mi, depending on the radius R and the cosmological constant L. This 
means that the form of the potential is given by Figs. [3^, and[3)D, respectively. Besides, the 
second derivative of the potential is always negative at m = m 2 . This means that the form 
of the potential is given by Fig. [3}o. 

Substituting equation (JBTl) into equation (J59l) we have 



d 2 V 



(R, m, L) 



1 



(R 6 L 2 + UR 12 - 20R 9 mL 2 + 2R 6 m 2 L 4 + 2L 4 mR 3 + 5L 4 ) . (59) 



dR 2 



2R 6 L 4 




(60) 
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FIG. 13: Case k\ = +1. The second derivative of the potential j^(R,m,L) calculated at m = mi 
(left) and at m = rri2 (right). We note that ^^(R,m = m\,L) can be positive or negative (the 
frontier between the two regions is given by equation pZ]) ) and that ^(R, m = m,2 , L) is always 
negative. 

and substituting equation (|58|) into equation fl59|) we get 

d 2 V,„ _ 3 



:(R,m = m 2 , L) = - or4Dfi x 



dR 2 v ' 

(-9i? 12 + 5R 6 V9R 12 + 2R e L 2 + 9L 4 + 9L 4 - L 2 V9R 12 + 2R e L 2 + 9L 4 ) . 

(61) 

Solving ^-(-R, m = m 1; L) = 0, we get 

L) « 1.816997838i? 6 . (62) 
Substituting equation ( 1571) into equation ( 1551) we have 

m = mi, L) = - 16 ^ R * x 
(-8L 4 R A + 18R 6 L 2 + 9R 12 + 3R 6 V9R 12 + 2R e L 2 + 9L 4 + 9L 4 + 3L 2 V9R 12 + 2i? 6 L 2 + 9L 4 ) , 

(63) 
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FIG. 14: Case k\ = +1. The potential V(R,m,L) calculated at m = mi (left) and at m = m,2 
(right). We note that V(R,m = mi,L) is always negative and that V(R,m = rri2,L) can be 
positive or negative. 



and substituting equation (1581) into equation (}55l) we get 
V(R, m = m 2 , L) = \ — - x 

(-8L 4 R A + 18R 6 L 2 + 9R 12 - 3R 6 V9R 12 + 2R 6 L 2 + 9L 4 + 9L 4 - 3LV9i? 12 + 2R e L 2 + 9L 4 ) . 

(64) 

We notice that V(R, m = m 1 , L) is always negative (see figure [Ml). Since V(R, m = m 2 , L) 
can be positive or negative, depending on the radius R and the cosmological constant L, we 
may have again formation of gravastars or black holes. 

VI. CONCLUSIONS 

In this paper, we have studied the problem of the stability of gravastars by constructing 



dynamical three-layer models of VW 3], which consists of an internal de Sitter spacetime, 
a dynamical infinitely thin shell of perfect fluid with the equation of state p = a, and an 
external Vaidya's spacetime. 

We have shown explicitly that the final output can be a black hole, an unstable gravastar, 
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R 



-0.1 - 



-0.2 - 



FIG. 15: Case k\ = +1. The potential V(R,m,L) calculated at m = mi, R c = 2 and L c = 5.75. 
This represents the formation of a "bounded excursion" gravastar. 

a stable gravastar or a "bounded excursion" gravastar, depending on the time evolution of 
the shell mass, the parameter L and the initial position R of the dynamical shell. 
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FIG. 16: Case k\ = +1. The potential V(R,m, L) calculated at m = mi, R c = 2 and L c = 
5.907256196. This represents the formation of a stable gravastar. 
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